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Abstract 

Identification by the Ulysses spacecraft of interstellar grains inside the planetary system provides 
a new window for the study of diffuse interstellar matter. Dust particles observed by Ulysses and 
confirmed by Galileo are more massive (> 10~^^ g) than the 'classical' interstellar grains. Even 
bigger grains (~ 10~^ g) were observed in form of interstellar meteors. We analyze the 
consequences of the plentiful existence of massive grains in the diffuse interstellar medium. 
Astronomically observed 'classical' interstellar grains can be described by a size distribution 
ranging from about 5 to 250 nm in radius (about 10~^^ to 10~^^ g). Lifetimes of these particles, 
due to mutual collisions in interstellar space, can be as short as 10^/ years, where / = 10 to 
1000, is the fraction of total lifetime to the time when grains are exposed to supernova shocks. 
Shattering is a source of the smallest of these grains, but grains more massive than about 10~^^ g 
of the classical interstellar grain population are rapidly destroyed. When applying the same 
shattering mechanism to the more massive grains found recently, we suggest that collisions of 
particles bigger than about 10~^^ g provide a source for smaller grains. Because massive grains 
couple to the interstellar gas only over large (100 to 1000 pc) length scales, the cosmic abundance 
ratio of gas-to-dust needs only to be preserved averaged over corresponding volumes of space. 
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1. Introduction 

Information about interstellar dust has been ex- 
tremely limited. Interstellar dust is usually astro- 
nomically observed by the obscuration (extinction) 
of starlight. Extinction of starlight, especially in the 
UV, provides a means to quantify the dust cross sec- 
tion along the line-of-sight between the observer and 
a distant star (e.g. Mathis, 1990). On a large scale, 
the starlight extinction is typically observable as av- 
eraged over galactic (kpc) distances. On a local scale, 
the amount of interstellar dust is indirectly inferred in 
the following way: the abundances of elements heav- 
ier than helium measured along the line-of-sight to a 
nearby star are compared with "cosmic abundances" 
and any missing heavy elements are assumed to be 
tied up in grains, which gives the amount of interstel- 
lar dust (cf. Fnsch, 1981; Fnsch et al., 1999). 

Galactic dust is believed to originate from a va- 
riety of different stars and stellar phenomena: e.g. 
carbon-rich stars, red giants, or supernovae, all of 
which provide dust with a different and characteris- 
tic chemical and isotopic signature (for a review, see 
Dorschner and Henning, 1995). A variety of preso- 
lar grains have been identified in primitive meteorites: 
e.g. diamonds, graphite, silicon carbide, or corundum 
grains {Zinner, 1998). However, the identified grains 
are only a minute fraction of the total material that 
went into the protoplanetary disk. The composition 
of the bulk of the grains is largely unknown. 

Interstellar dust grains passing through the plan- 
etary system have been detected by the dust de- 
tector onboard the Ulysses spacecraft [Griin et al., 

1993) . These observations provided imique identifi- 
cation of interstellar grains by three characteristics: 
1. At Jupiter's distance, the grains have been found 
to move on retrograde trajectories, opposite to the 
orbits of most interplanetary grains, and this flow di- 
rection coincides with that of interstellar gas ( Witte 
et al., 1993), 2. a constant flux has been observed 
at all latitudes above the ecliptic plane [Griin et al., 
1997; Kriiger et al., 1999), while interplanetary dust 
displayed a strong concentration towards the eclip- 
tic plane, and 3. the measured speeds (despite their 
substantial uncertainties) of interstellar grains were 
higher than the locel Solar System escape speed, even 
if one neglects radiation pressure effects {Griin et al., 

1994) . 

It is known that neutral interstellar gas flows through 

the solar system with a speed of 26 km from 
the direction of 253° ecliptic longitude and 5° ecliptic 



latitude {Lallement, 1993, Witte et al., 1993). The 
flow of interstellar dust detected by Ulysses coin- 
cided with this direction [Frisch et al., 1999), and 
persisted during Ulysses' tour over the poles of the 
sun. Interstellar dust was identified as close as 1.8 
AU from the sun at ecliptic latitudes above 50° ( Griin 
et al, 1997). In addition, dust impact measurements 
from the Hiten satellite, in high eccentric orbit about 
the Earth, gave indications that interstellar mete- 
oroids indeed reached the Earth's orbit {Svedhem et 
al, 1996). Other dust measurements by Galileo con- 
firmed these findings {Baguhl et al., 1995; Baguhl et 
al., 1996). Modeling of the Galileo and Ulysses data 
suggested that up to 30% of dust flux, with masses 
above 10~^^ g, at 1 AU is of interstellar origin [Griin 
et al, 1997; Landgraf, 1999). 

The masses of clearly identified interstellar grains 
range from 10^^^ g to above 10~^^ g with a maxi- 
mum at about 10~^^ g. More recently, even bigger 
(» 10~^ g) interstellar meteors have been reliably 
identified [Taylor et al, 1996, Baggaley, 1999). The 
deficiency of small grain masses (< 10~^^ g) is not 
solely introduced by the detection threshold of the 
Galileo and Ulysses instruments, but indicates a de- 
pletion of small interstellar grains in the heliosphere. 
Estimates of the filtering of 0.1 /im-sized and smaller 
electrically charged grains in the heliospheric bow 
shock region and in the heliosphere itself [Frisch et 
at, 1999; Linde and Gombosi, 1999; Landgraf, 1999), 
show that these particles are strongly impeded from 
entering the inner planetary system by the interac- 
tion with the ambient magnetic field that is carried 
by the solar wind. We will assume in the following 
analysis that, in addition to the big particles observed 
within the heliosphere, there exist also the astronom- 
ically observed smaller grains in the local interstellar 
medium. 

Astronomically observed interstellar dust is con- 
veniently described by the MRN size distribution 
[Mathis, et al, 1977) ranging from 5 to 250 nm in 
size (about 10~^^ to 10"^'^ g). The size distribution 
of these particles can be described by a power law 
distribution (with an exponent of -3.5) and the to- 
tal amount is related to the volatile density in the 
local cloud. On the basis of this size distribution, 
Jones et al. (1996) find that collisional shattering is 
a source of the smallest of these grains, but grains 
bigger than about 30 nm are rapidly depleted, much 
faster than they can be replenished by condensation 
in stellar outflows [Jones et at, 1994). In this study 
we investigate whether collisions with the newly found 



3 



bigger particles could provide a source for grains big- 
ger than those generated by collisions from the classic 
interstellar grains alone (cf. Jones et al., 1997). 

The mass distribution of interstellar grains mea- 
sured by Galileo and Ulysses only overlaps with the 
biggest masses of the 'classical' MRN distribution, 
and extends to much bigger particles. For this analy- 
sis, we will assume that the mass distribution of the 
big particles extends to the mass of the observed in- 
terstellar radar meteors (« 10~^ g, see Landgraf et 
al, 1999, for a discussion of the mass distribution ob- 
served by spacecraft and radar observations). While 
the grain flux up to 10"^'^ g is well characterized 
by spacecraft measurements the meteor flux at about 
3 X lO"*" g has still large uncertainties. Landgraf et 
al. (1999) give a cumulative mass flux with a slope 
of —1.1 ± 0.1 (corresponding to a —4.3 ± 0.3 slope 
of the size distribution) for the Ulysses and Galileo 
measurements alone, however, the radar data suggest 
a significantly flatter slope of the combined mass dis- 
tribution, and a slope of -4 is compatible with both 
measurements. These big particles are mostly hidden 
from astronomical observations because their contri- 
bution to the total cross section for extinction is low. 

Where do big grains come from? Since big grains 
are difficult to observe, here we can only speculate. 
Stellar outflows (star dust) and explosion shells dis- 
play only small grains that are injected into the inter- 
stellar medium. Circumstantial evidence for micron 
sized and bigger grains outside of our own solar sys- 
tem exists mostly in dust disks around Vega type stars 
(cf. eg. Backman and Paresc, 1993). It is in molec- 
ular clouds, star forming regions and, especially, in 
proto-planetary disks where one would expect grain 
growth, however, the mechanism by which these large 
grains escape into the diffuse medium remains un- 
clear. Excess emission at millimeter wavelengths from 
very cold dust in the Galactic plane is interpreted by 
Rowan-Robinson (1992) to be due to particles of sev- 
eral 10 /Ltm in size. 

In the main part of this paper, we will demonstrate 
that the mere existence of a significant number of big 
particles in parts of the diffuse interstellar medium 
has locally profound consequences for the evolution of 
interstellar material. The big particles' masses pro- 
vide a significant collisional reservoir for smaller par- 
ticles, in those interstellar regions where big particles 
are abundant. However, it is not the scope of this pa- 
per to develop an evolutionary model of interstellar 
dust in the local interstellar medium, especially, since 
source terms for interstellar grains are not considered. 



We will make simplifying assumptions for the size 
distribution and the collision speed in order to arrive 
at an estimate of the collisional effects of big inter- 
stellar grains. First, we assume a size distribution 
of big grains that is compatible with that observed 
by Ulysses, and combine it with a standard model of 
locally unobserved small classic grains. The conse- 
quence of this assumption will be that the results ob- 
tained are strictly valid only in regions of space where 
such a grain size distribution is found. Furthermore, 
we will assume a constant speed of 100 km s^-'^, as 
the effective collision speed at which particles of dif- 
ferent sizes collide. This speed is derived from the 
average speed that supernova shocks travel through 
the interstellar medimn (cf. Jones et al., 1996). The 
collisional effects arc mitigated by the factor f that is 
the fraction of total lifetime when grains are exposed 
to shocks, where / = 10 to 1000. While the abso- 
lute numbers may be uncertain by a factor 10, the 
qualitative results that we emphasize are still valid. 

The mass added by the big grains to the local in- 
terstellar medium is not violating well founded cosmic 
abundance arguments. In the discussions portion of 
the paper, we will argue that the dust-to-gas mass ra- 
tio may vary locally, and that the cosmic abundance 
needs only to be preserved from averages over large 
volumes of space. 

2. Size and mass distributions of 
grains in the local interstellar medium 

The MRN size distribution of grains in the dif- 
fuse interstellar medium is given by a power law 
dn = AnHa~°'da, with slope a = 3.5, and A = 7.76 x 
lO"^*" cm~2-^ per H nucleus (Mathis et al, 1977). In 
the following, we will assume uh = 0.3 cm~^, a typ- 
ical value for the local interstellar medium (Prisch et 
al., 1999). The grain sizes, a, range from 5 to 250 nm 
(m = 1.7 X 10-18 g to 2.2 x lO'^^ q^j. 
cussion, we follow a similar discussion by Griin et al. 
(1985) for the distribution of interplanetary dust but 
we use, instead of the size distribution, the mass dis- 
tribution after assuming spherical particles of density 
p = 3 g cm"^. For mass distributions covering a wide 
mass range, it is convenient to use the logarithmic 
differential distribution: 

dn = CMKNm^~ dlogm (1) 

with the constant 

Cmrn = Ann ( ) ^InlO. 



4 



log. o (cm) 




-18-16-14-12-10-8 -6 
log. m (g) 



■18-16-14-12-10-8 -6 
log. m (g) 



■18-16-14-12-10-8 -6 
log. m (g) 



Figure 1. Mass distribution and moments of the mass distribution (logarithmic differential number-, cross 
section-, and mass-density distributions) for the MRN (bold dashed) and the extended distributions. 



The new extended mass distribution consists of 
three parts: a MRN-type distribution (slope —a = 
—3.5) for small particles (about 10^^^ g to 10"^'* g), 
a distribution with a steeper slope {—a — —4.0) for 
bigger particles (about 10^^^ g to 10~^ g), and a tran- 
sition region in between. We have assumed that the 
slope of the big particle distribution holds up to the 
size of the radar meteors. The extended mass distri- 
bution is represented by 
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where —(iS ~ —3.5 is the slope 

for small particles and — /3 = —4 is the slope for big 
particles, and 5 — 0.5 controls the transition. Again, 
we transform this extended distribution to the loga- 
rithmic differential distribution: 
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the transition mass mt = 1.5 x 



g and Ccxt = 9.0 x 10" 



The extended 



distribution is valid from m = 1.7 x 10 g to 2.2 x 
10^^ g (a = 5 nm to a = 25 /im). 

The cross-section and the mass distributions are 
given by 
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Fig. 1 shows the MRN and extended mass dis- 
tributions and their moments: cross section den- 
sity ua and mass density distributions. For 
the extended distribution, the total mass density is 
M = J Umdlogm = 1.78 x lO'^^ g cm'^ (MRN: 
4.97 X lO"^*" g cm~"^) and the total cross section den- 
sity is ^ = J HAdlogm = 3.90 x 10~^^ cm^ cm"'' 
(MRN: 3.63 x 10"^^ cm^ cm'^). The cross section 
density of the extended distribution is not signifi- 
cantly increased (-1-7%) over that of the MRN distri- 
bution and, consequently, the interpretation of extinc- 
tion measurements at visible and UV wavelengths is 
not affected. However, most mass density is contained 
in the big particles (the total mass of the extended 
distribution is about a factor 3.5 increased over that 
in the MRN distribution). 
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3. Collision Dynamics 

We now determine the effects of mutual collisions 
in a cloud of interstellar matter like that currently 
surrounding the solar system. On an average of ev- 
ery few 10^ years, a supernova shock passes through 
the diffuse medium with relative speeds of about 50 
to 200 km s~^ {Jones et ai, 1996). The pressure 
jump and the entrained magnetic field first acceler- 
ates the gas then betatron acceleration causes the 
small grains, and eventually the bigger grains, to 
reach their post shock speed. In this way, relative 
speeds between grains of various sizes are introduced. 
In this study, we will focus on the coUisional de- 
struction and the generation of fragments by colli- 
sions between grains, and we will ignore sputtering of 
grains by the fast moving gas, which may addition- 
ally erode small grains in very fast shocks. We assume 
a quasi-stationary situation, by taking into account, 
only the shock environment where the relative speeds 
between the grains are of the order of the shock speed. 
Hence, we will not follow the coUisional time evolu- 
tion of the size distribution, and therefore we assume 
that the grain size distribution does not change with 
time. Taking into account that interstellar matter cy- 
cles repeatedly through a warm and cold interstellar 
medium (where it is effectively shielded from shocks), 
the effective lifetimes may be a factor / = 10 to 1000 
times longer than calculated. 

In the subsequent section, we follow the collision 
formulation of Griin et al. (1985) for interplanetary 
dust and that of Jones et al. (1996) for interstellar 
dust. An important parameter describing the effect 
of collisions is F, the ratio of the target and projectile 
masses, at which a "catastrophic" collision occurs, i.e. 
a collision that shatters the target particle completely. 
Jones et al. (1996) derives this value T = 10^ to 10'') 
from first principles for interstellar particles in a shock 
environment. Griin et al. (1985) used values F = 10^ 
to 10^) for interplanetary collisions that were derived 
from laboratory experiments with cm-sized projectiles 
and basalt targets at impact speeds below 10 km s""'. 
Since the Jones et al. values were derived for the size 
and speed range applicable to interstellar grains, we 
use their value (F = 3 x 10'^). 

During the passage of a supernova shock through 
the interstellar medium, the relative speeds between 
particles of different sizes vary over a wide range with 
the maximum relative speed close to the shock speed 
itself. Any target particle of a given mass is most 
effectively destroyed by projectiles that are up to a 



factor F smaller than the target particle. Relative 
speeds between different massive particles reach shock 
speed. Even particles of the same mass can get high 
relative speeds, since shock acceleration is facilitated 
by the pick-up process, where charged grains gyrate 
about the moving magnetic field and hence colliding 
particles can be at different phases of their gyration. 
CoUisional effects are strongly dependent on the col- 
lision speed (see e.g. Jones et al, 1996), and the 
highest speeds have the strongest effect. Therefore, 
we assume v — 100 kms^^ as the effective speed of all 
catastrophic collisions during the passage of a shock. 

The collision rate c(mi) of target particles of mass 
ml is given by 

c(m\) = - / <j{mi,m2)v— alogm2 

Jmi/r dlogm2 

(6) 

where a {mi, 1712) is the coUisional cross section. 
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Figure 2. CoUisional lifetimes of grains in the lo- 
cal interstellar medium (nn = 0.3 cm.~^) for projec- 
tiles from the MRN (bold dashed) and the extended 
distributions. 

From the collision rate we can determine the mean 
coUisional lifetime T{m) of grains of mass m: 

r{m) = (7) 
c(m) 

Fig. 2 shows that coUisional lifetimes for the 
biggest particles of the MRN distribution (10"^'^ g) 
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are of order a few 10^ years, giving rise to a fast de- 
pletion of these particles, and a steepening slope of 
the MRN distribution as described by Jones et al. 
(1996). The biggest particles of the extended distri- 
bution (10~^ g) have lifetimes about 100 times longer 
than those of the biggest MRN-particles. 

In the next step, we calculate the mass density 
destruction rate dMd shattered in catastrophic col- 
lisions: 



dMd{mi) 
dlogmi 



dn{mi)mi 



dlogmi 
dn{m2) 

dlog TO2 



(T(mi, m2) 



i/r 

dlog 7712 . 



(8) 



It must be noted that the mass of the colliding 
particles is not lost, but reappears in the distribution 
of fragment particles. From laboratory experiments, 
the fragment mass distribution is given by 



dG{m, nil, 1712) 
dlog m 



— const. • m ^ , for m < ttil 



(9) 



with rj — 0.77 (after Jones et ai, 1996), rriL is 
the mass of the largest fragment (we will use a value 
of 5 X 10"'^ of the target mass) and the constant in 
equation (^ can be evaluated from the conservation of 
mass (assuming no significant losses by evaporation). 

The mass density generation rate Mg of fragments 
can be calculated from 
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Figure 3. Mass distributions of destroyed parti- 
cles and the generated fragments for the MRN (bold 
dashed) and extended distributions. 

Collisions are an effective source of small grains 
(cf. Jones et al, 1996): more fragments are gener- 
ated than grains are destroyed in each mass interval 
below 10" 16 g and 10" g for the MRN and the ex- 
tended mass distribution, respectively. The situation 
is reversed above these thresholds, however, while big- 
ger grains are rapidly lost from the MRN distribution, 
the losses from the extended mass distribution are re- 
duced by coUisional fragments that are generated up 
to 10-1° s. 
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(10) 



where /i = 10 1* g is the mass of the smallest particle 
considered as projectiles. 

Fig. 3 shows the mass distributions of the de- 
stroyed particles and the generated fragments for the 
extended distributions. The total mass shattered 
(equals the total mass of generated fragments) per 
unit time and volume for the MRN distribution is 
5.7 X lO-** rmg m"^ s-^) whereas the processed mass 
for the extended distribution is 1.1 x 10"'^ g m"'^ s^^). 



4. Discussion 

A second consequence of the existence of big in- 
terstellar grains in the local interstellar medium is 
that the gas-to-dust ratio may deviate significantly 
from that derived from cosmic abundance arguments 
(about 1% of the total mass is in refractory dust; cf. 
Frisch et ai, 1999). In this section, we show that 
massive dust grains couple to the interstellar gas over 
much larger distance scales than do the small classic 
interstellar grains. As a consequence, the dust-to-gas 
mass ratio may vary locally (on scales of a few pc), de- 
pending on the relative abundance of massive grains, 
and cosmic abundance arguments are valid only av- 
eraged over large volumes of space (100 to 1000 pc 
dimension). In the following section, we calculate rel- 
evant length scales that govern grains in the diffuse 
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interstellar medium and discuss consequences for the 
gas-to-dust mass ratio. 

For comparison with other processes that affect 
grains in the diffuse interstellar medium, we calculate 
the collisional length scale: 



coll 



(11) 



Again we assume that dust particles move with 
a speed of u = 100 kms~^ with respect to the am- 
bient medium. This length scale is a lower limit - 
it could be a factor f longer. Next we calculate the 
length scales over which interstellar dust couples to 
other constituents of the ambient diffuse interstellar 
medium: gas and fields. The frictional scale Idrag 
over which the dust couples to the gas (cf. Morfill 
and Griin, 1979) is given by the path length it takes 
for the grain to sweep-up its own mass, md, in form 
of interstellar gas (hydrogen, rriH = 1.67 x lO^^'' g): 



Id 



rag 



AdUHmH 



(12) 



An important scale for an interstellar grain is its 
coupling to the ambient magnetic field. We will 
assume that interstellar particles in the local inter- 
stellar medium are charged typically to U = Q.5V 
{Griin and Svestka, 1996). Therefore, their charge 
is given hy q = AweoUa, with permittivity, eo = 
8.859 X 10"^^ CV^^ m"^, and a in meters. The gyro 
radius, Igyr, in the ambient magnetic field, B = 0.5 nT 
(estimated average value in the local diffuse interstel- 
lar medium, Holzer, 1989), is 



mv 



(13) 



Fig. 4 compares the scales for collisional de- 
struction, gas friction, and gyration, assuming that 
the local condition are representative for the dif- 
fuse interstellar medium over a much larger volume. 
The shortest scale is the interaction with the mag- 
netic field. If we further assume that the interstel- 
lar magnetic field is coupled to the gas via the ion- 
ized component, the gyro radius is the shortest scale 
over which the dust couples to the gas. Neverthe- 
less, the biggest particles (10~^ g) can travel over 
about one kpc before they are captured by the mag- 
netic field, and therefore, big particles may be gen- 
erated several 100 pc from the place where they are 
found. It is evident that big particles couple over 
much larger scales to the ambient diffuse interstel- 
lar medium (gas and fields) than small particles do. 



therefore, they may be unrelated to the local gas 
density and the small particles except via collisions. 
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Figure 4. Length scales (1 pc — 3 x 10^^ m) for 
collisional destruction (extended mass distribution 
only), gas friction (gas density Uh = 0.3 cm~^), 
and gyration (surface potential U = 0.5 V, magnetic 
field B = 5nT). 

Consequently the dust-to-gas mass ratio may vary 
locally, and hence, the cosmic abundance (about 
lonly over large scales (100 pc to 1 kpc). Another 
consequence of the much longer coupling length of 
big particles is that there is everywhere a significant 
relative speed (of the order of 30 kms^^) between big 
grains and smaller grains that couple to the local gas. 
Initial estimates indicate that the collisional effects of 
this constant grinding may be as important for the 
generation of small particles as those induced by su- 
pernova shocks. 

The most important consequences of the big parti- 
cle population in the local diffuse interstellar medium 
are that: 1. most mass is in big particles, 2. massive 
particles have long collisional lifetimes ( 10^ years), 
3. big particles provide a source for small MRN-type 
particles as long as big grains are present, 4. massive 
grains couple to the gas over length scales of 100 to 
1000 pc, 5. the dust-to-gas mass ratio varies locally 
depending on the relative amount of massive parti- 
cles, 6. cosmic abundance ratio of gas to dust is only 
valid over kpc distances. 

The purpose of the paper is to draw the attention 
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to effects of big interstellar grains that have not been 
previously recognized. A full evolutionary model of 
interstellar dust in the diffuse interstellar medium is 
left to future work. 
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